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SUMMARY 

Transport equations for a partially ionized gas in an electric field are 
derived from the Boltzmann equation by the Grad 13 -moment method. Charged 
particle interactions are described by Coulomb forces, vhile all other particle 
interactions are described by a rigid elastic sphere model. First approxima- 
tions to currents, stresses, and heat fluxes are obtained. The resulting 
transport relations are suitable for use in situations vhere the electrons have 
an elevated temperature. 


INTRODUCTION 

There are two methods currently in use in the description of transport 
phenomena in multicomponent gas mixtures. The Chapman-Enskog method (refs. 1 
and 2) produces a set of approximations to currents, stresses, and heat flovs. 
The theory is restricted, however, to small departures from local equilibrium; 
in particular there is no provision for assigning separate temperatures to the 
components of the mixture. The Grad method (refs. 3 to 8) results in a set of 
partial differential equations for the coefficients in the expansions of the 
species distribution functions about local Maxwellians. The first few expan- 
sion coefficients are simply related to currents, stresses, and heat fluxes. 
This theory permits greater departures from equilibrium than does the Chapman- 
Enskog solution and, moreover, permits the assignment of different temperatures 
to each component. For small departures from the local Maxwellians, a series 
of approximations can be generated for currents, stresses, and heat flows from 
the differential equations. 

In the presence of an electric field, the electrons in an ionized gas are 
usually characterized by a higher temperature than the heavier particles. Pre- 
vious treatments have been limited to either such low degrees of ionization 
that only electron- neutral and ion-neutral interactions need be considered 
(refs. 9 and 10 ) or to such high degrees of ionization that the gas may be 
treated as fully ionized (refs. 8 and ll). In the present work, the descrip- 
tion of the general n-component mixture obtained by the Grad method (refs. 4 
to 7) is worked out in detail for a three- component mixture of electrons, ions, 
and neutral atoms. A first approximation to currents, stresses, and heat flows 
is obtained through an iterative scheme. 



There is some question as to the applicability of the standard collision 
integral (see eq. (2)) in the case of Coulomb interactions because of the long 
range nature of these forces. Comparisons with the usually preferred Fokker- 
Planck treatment, however, show good agreement. The second approximation to 
the electrical conductivity as obtained by the Chapman-Enskog method differs 
by only 2.1 percent from the Fokker-Planck value (refs. 11 and 12), while the 
value obtained by the 13-moment method (ref. 8) differs by 2.6 percent from 
that obtained with the Fokker-Planck equation. The problem of collective phe- 
nomena, such as the interaction of a high-speed electron with its wake of 
plasma oscillation (ref. 13), is not included in the previous treatments and 
has been neglected herein. 

Cartesian tensor notation is used throughout this report. Greek sub- 
scripts are used for tensor indices with summation on double subscripts. 
English superscripts refer to species. 


THEORY 


The theory has been treated in detail elsewhere (refs. 3 to 7). A brief 
resume is given here. The distribution function f^ for each species is 
assumed to satisfy the Boltzmann equation 
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where is the particle velocity, the position, and E a the electric 

field. The collision integrals 1*3^ are given by 

I Jk = y [f J (T , )f k (T{) - f j (T)f k (?!)]gb db de djfi (2) 

where g = |q - | , b is the impact parameter, e the azimuthal angle, and 
-> -» -» -» 

| are the final values of the velocities Ijiq after an encounter. (All 
symbols are defined in appendix A. ) 

The are expanded in Hermite polynomials (ref. 14) about local 

Maxwellian distributions fJ (*-*). In the 13-moment approximation, f «3 takes 
the form 


fj(°) n 

L + a^ 1 )v a + \ a^ 2) v a v p + ^ aJ (3) v a (v 2 - 5)] 

(3) 

v a ^ 

(a.j) / c a a J - c a - - w a 
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hj = 4 

(5) 
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( 6 ) 


fj(0) „ exp(.| v2) 

a^*^ = (aJ) 1 / 2 (7) 

a a 



In these equations , u^ is the drift velocity of species j with respect to 

CL • ^ 

the mean mass velocity w a , p^p is the nonhydrostatic stress, and is a 

reduced heat flux. 

In order for the approximation of equation (3) to be valid, the expansion 
coefficients must be small compared to unity. In particular this requires that 

« (kTJ/mJ)' 1 -/^ and pJ « p^ 

CG CGp 

Equations for n J , p J , u^, p^p, and are generated by taking moments 

of the Boltzmann equation. The higher moments and the collision integrals 
occurring in these equations are evaluated in terms of preceding quantities by 
means of the approximation to fJ in equation ( 3 ). The result for an 
m-component mixture is 13m ^ coupled ^ partial differential equations in the 
13m dependent variables n^ , pJ, uJ, P^* and - 

The collision integrals have been considered in detail in the literature 
(refs. 3 to 7), and the results for elastic collisions are listed in convenient 
form in reference 5. For the Coulomb interaction between charged particles, 
the integration over the impact parameter b has been cut off at the Debye 
length (ref. 15): 



The integration over velocity is taken from reference 5 with a small correction 
(see appendix B). The rigid elastic sphere model has been used for all other 
interactions. The collision terms, in general, are exceedingly complex. The 
following restrictions, however, result in a considerable reduction in com- 
plexity for a plasma: 
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m e « m a m 1 = m a 

Furthermore, as in references 5, 7, and 8, all terms quadratic in u^, pj^, RJ 
have been neglected. 


MOMENT EQUATIONS 

. . * 2 

When the moment equations for nr 1 , m°c a , and l/2 m J c are summed over 
the species, the conservation equations for the mixture are obtained: 


Continuity: 
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Momentum: 
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Energy : 
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where 
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The 13 moment equations for each species are generated from m J , m J Cy., 
\ m^c 2 , and m^ ^yc^ - | c 2 5 r ^ ^ m ^ G y |c 2 - —j 

Continuity: 
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( 19 ) 


£§£ + n J ^ P 

Dt dx a 


Momentum: 


mJ| 


DP: 


Dt + r r 3x“ 1 a ^ 
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= - ^ Tp^n k rj - n^r k ^ + ^ a^ k ^n k Rj - ^ nJR k ^ (20) 


k k 

where equation (13) is used to define a new force 
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Energy : 


Stress : 


^ + „J + i Ml . 2 

Dt dx^, 5 tef 




+ 2<P^ 


. dw-> 


“rS"( 


= - D 


where for any second-rank tensor 

{Axp} s | (Axp + - | A rr S ap') 


Heat flux: 


S Pa 


’ar .5 wpJ V& 
~ 1 + " ^ 


(21) 


m (I r J ) + I pJ ST + ^ ■ F « r » + P aS ^ = - E S J Vn k (TJ - T“) (22) 


(23) 


(24) 


-i dwx ■? dwy- k (l i 5 -i i - raT . -. 

dT + ®r + K a ^ + ^ (-2 ^ ^ + 2 pJ ^ + TJ + 2 “V 

- ^ * *»♦! {tk + *Z k +q ± -£ v3k -f < 25 > 


Collision coefficients for the system of electrons, ions, and neutral atoms 
are listed in appendix B. 
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In writing the collision terms as they appear on the right sides of equa- 
tions (23) and (24) and in appendix B, an order- of -magnitude analysis has been 
used to discard several terms. The complete analysis is too complicated to be 
reproduced here, but a sample portion is given in appendix C. The end result 
is that the collision terms appear in their given, forms, and the electron 
stress and heat flux equations uncouple from those for the ions and atoms pro- 
vided that certain inequalities are satisfied. The most stringent restriction 
is given by 


z Off; /Tf ^ 1 / 2 /^ 1 / 2 

Q ia VW Vm a / 


« 1 


(26) 


The collision coefficients corresponding to the neglected terms are set equal 
to zero in appendix B. 


TRANSPORT APPROXIMATIONS 


It has been proposed that approximations to currents, stresses, and heat 
fluxes be generated from the moment equations by an iterative scheme (refs. 1, 
6, 7, and 8). Let T and L be a macroscopic time and length characterizing 
a flow. It is assumed that T“1 and wL”l are small compared with the re- 
laxation frequencies appearing as coefficients on the right sides of the moment 
equations. In addition, it is assumed that gradients of currents, stresses, 
and heat flows are small. Then, . as a first approximation, a set of algebraic 
equations is obtained for P.j, pj^ and R^ as functions of nJ, tJ, Ey, and 

Wy. and their gradients. Higher approximations are obtained by substituting 
the value derived in the previous approximation into those terms originally 
neglected on the left side of the equations. 


The first approximations are now considered in detail. Under the preced- 
ing assumptions equations (20), (23), and (25) become 


- + n jpJ = 


dXy 


ipj _ l _ n ^ qjk |nk 4 - i -A (27 > 


+ 2 P*4 - E ffiJk & 


(28) 


- I pJ i i; + fi E °? 3 r°- - E - E r ^( nkr r - nJr ?) (29) 


with 
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' p OXy 


(30) 
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The vector equations (27) are not independent. As the PJ. are referred to the 
mass velocity Wy 

(31) 


s 


m k r k = o 


Inverting equations (28) gives 


- 2 ” e ( e ¥t} 


v- - - 21,1 Wr + 2,1 


V - " 2 ^M - ^ 


i {zeE r r5;} 
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(52) 

(35) 

(34) 


where the viscosities \± are given by 
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(35) 


and the coefficients by 



Jt e = (<D ee )”^ 





1 

Tt^~ = 03 aa (a3^^a3 aa - 03^ a 03 a ^) 

> 

(36) 


Jt a = “fi Jt 1 
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Inverting 

equations (29) gives 




-U 

II 

>> 

CD 
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t) + 

D ea (n a r a - n e r a ) 

(37) 


^ = - Al S + | TlE a4a + Bi ' 

(n^ 

- « ±r r) 

(38) 
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(39) 


*? ’ - AS If - § T X.4a + Da (“ ll r - 

where the thermal conductivities A^ are given by 

fs ^ 

_ 5 _k_ p 
2 m e v ee 

_ 5 _k_ v 8 ^) 1 - v ia p a ^ 

2 j^a yiiy aa _ v 

A a = — — v^p a - v a ^p^ 

^ m a v i:L v aa - v^ a v a ^ 


the coefficients by 
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and the coefficients D by 
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v xx y aa 
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(40) 
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(42) 


The final step is the inversion of equations (27) with equation (31) used 
in place of the ion equation. Equations (37) to (39) are used to eliminate the 
reduced heat fluxes R^. The following diffusion equations result: 




(43) 
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where the forces dy are defined b y 

_e 

-T 


n e d® = - n e F® + (a ei ni + a.ea n a) J T e E a pJ a ) 




dxy 


(n 1 + n a )dS = 4_ _ n a F a - a ea n a (A e 4— + -% T e E p® 4 
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and the transport coefficients are given by 
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It has been assumed that 


n e « n j 
m a 


m ia ^ ea 


m‘ 
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Tp 0 ' « rj ea n e T} ea « n 1 !) 13. 


(48) 


In general the term in d a in equation (43) is negligible. 


SUMMARY OF RESULTS 

The 13 moment transport equations have been obtained for a partially 
ionized gas in an electric field. All particle interactions have been 
described in terms of binary collisions; the shielded Coulomb model has been 
used for those between charged particles and the rigid elastic sphere model 
for all interactions involving neutral atoms. According to Grad these equa- 
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tions should be valid even in cases -where there is an appreciable variation in 
properties over the distance of a mean free path. 

First approximations to currents , stresses,, and heat fluxes have been ob- 
tained from the 13 moment equations. These approximations include expressions 
for diffusion coefficients, viscosities, and thermal conductivities and exhibit 
cross -coupling between currents and heat fluxes. Hie resultant transport rela- 
tions are suitable for use in situations where the electrons have an elevated 
temperature. 

National Aeronautics and Space Administration 
Lewis Research Center 

Cleveland, Ohio, February 19, 1964 
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APPENDIX A 


a 1 

a 


j(n) 


CO- 

D ' 5 
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E, 


’a 
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f j 
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h 

iJk 

Ja 
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n J 


pJ 

ap 

pj 


*ap 

Q jk 


SYMBOLS 

expansion coefficient in distribution function, eq. (3), where 
n = 0, 1, 2, 3 

impact parameter 

peculiar velocity, g a - w a 

coefficient, eqs. (42) 

diffusion force, eqs. (46) 

electric field 

electronic charge 

charge on particle of species j 

force, eq. (21) 

distribution function of species j 

magnitude of relative velocity 

Debye length, eq. (10) 

collision integral, eq. (2) 

electric current 

Boltzmann constant 

mass of particle of species j 

number density 

stress tensor 

pressure, l/3 

nonhydrostatic stress tensor, P^p - p J 6 a p 
rigid sphere collision cross section 
heat flux 

reduced heat flux, q^ - 5/2 
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rpJ 

t 


U 


a 


v a 


w a 

x a 

Z 

a jk 

pjk 

r J 

a 

r <jk 

Sap 


n jk 

e«J 

fcJ k 

A 

aJ 

p.j 

vJ k 

it j 

P 

a 

tJ 


kinetic temperature of species j 
time 

diffusion velocity 

dimensionless peculiar velocity 

mean mass velocity 

position vector 

ion charge in units of e 

coefficient, eq. (20) 

coefficient, eq. (22) 

diffusion current, n^u^ 

oo 

coefficient, eq. (25) 

Kronecker delta 
azimuthal angle 
coefficient, eqs. (47) 
coefficient, eq. (20) 
coefficient, eqs. (47) 

collision parameter, eqs. (B22) to (B25) 

parameter, eqs. (B26) to (B28) 

thermal conductivity, eqs. (40) 

viscosity, eqs. (35) 

coefficient, eq. (25) 

particle velocity 

coefficient, eqs. (36) 

mass density of mixture 

charge density of mixture 

coefficient, eqs. (4l) 
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coefficient, eq. (23) 
Subscripts : 

tensor indices 
Superscripts : 
a atom 

e electron 

i ion 

any species 

( 1 ) velocities after collision 



APPENDIX B 


COLLISION COEFFICIENTS 


The collision coefficients appearing in the species moment equations for 
the electron, ion, neutral-atom system are as follows: 


T) ei - = Ti le = i y/2. Z 2 m e K ee (ln A - r) 
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co aa = 2 ± K ia n i + Bb K aa n a 
15 5 


84: _j_ 16 

(cD ei ~ 0^- e ~ CD ea ~ CD' 


ae 


0) 


(B13) 


14 


I 



v ee = K ee |i§ /in A ee - x + — )n e + -£■§ Z 2 /ln A - y — — V 4 ! + — K ea n a (B14) 

115 \ 2/ 15 \ 13/ J 15 

v 11 = M z 4 ^ l/2 ^ 3/2 *ee(i n A ii _ r + !) n i + || K i a n a (B15) 

via = _ | K ia n i (B16) 
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v ai = - 1 K ia n a ( B17 ) 
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v aa _ 59 K ia n i + || /<aa n a (B18) 


( v ei « v ie v ea « v ae « q) 


T e ^ = 2 ^2 Z^kT e K ee ^ln A - y - 
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Y*ea _ _ 4 
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A _ 2hkT e 
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A ee = ZA 
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(B28) 


Ai 1 . JEl 

ZT e 

and x is Euler's constant,, 0.577. 

It should he noted that equation (A6) of reference 5 should he corrected 
to read 




With this correction the last two equations in section (5-3) of refer- 
ence 5 become 
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APPENDIX C 


SIMPLIFICATION OF HEAT-FLUX EQUATIONS 


In the evaluation of the collision integrals in the ion heat-flux equa- 
tion these two terms appear: 

f V 5 z2 $ A - r)n<^ 


and 


VS/.evVs 


15 


16 K eef 

W7 \T a / l 


K ee ( In - r + 


i ) ni 4 


(Cl) 


The first term will be negligible compared with the second if the following 
inequality is satisfied: 


« 1 (C2) 


\l/2/ T a\3/2 

p e u^ 

7 W 

4 


where p e = n e kT® and In has been assumed to be of the same order as 

In A. 


For the fully ionized gas (n^ » n a ), the ion current is very small, and 
the inequality is assumed to hold (ref. 5). For nl ~ n a and nl « n a , 
equations (38), (42), (45), (B15) to (B18), (B2l), (B22), and (B24) give 


P e ul 

T 


p e y i:L v aa „ y ia y ai kT a 

pi v aa + v ia r ia (n a + n 1 ) 

pf v ii kT a 

pi yia( n a + n i) 



Under conditions in which the moment equations are assumed to hold n e « Zn^. 
Equations (C2) and (C3) lead to these two inequalities in the two density- 
ranges: 
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« 1 


(C4) 



(C5) 


These relations have heen presumed, to hold. 

Analyses similar to this one have "been used to simplify the other moment 
equations. 
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